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Abstract 

We study the handbag contribution to two-photon annihilation into baryon-antiba- 
ryon pairs at large energy and momentum transfer. We derive factorization of the pro- 
cess amplitude into a hard 77 — * qq subprocess and form factors describing the soft 
qq — * BB transition, assuming that the process is dominated by configurations where 
the (anti)quark approximately carries the full momentum of the (anti)baryon. The form 
factors represent moments of time-like generalized parton distributions, so-called BB dis- 
tribution amplitudes. A characteristic feature of the handbag mechanism is the absence 
of isospin-two components in the final state, which in combination with flavor symmetry 
provides relations among the form factors for the members of the lowest-lying baryon 
octet. Assuming dominance of the handbag contribution, we can describe current experi- 
mental data with form factors of plausible size, and predict the cross sections of presently 
unmeasured BB channels. 
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1 Introduction 



In this article we study the annihilation of two photons into baryon-antibaryon (BB) pairs at 
large Mandelstam variables s ~ — t ~ — u in the handbag approach recently developed for two- 
photon annihilation into pairs of mesons 0. As in the meson case, the handbag amplitude (see 
Fig. 0) factorizes into a hard 77 — > qq subprocess and form factors representing moments of 
generalized distribution amplitudes |2|, fj . These distribution amplitudes are time-like versions 
of generalized parton distributions, which encode the soft physics information in processes such 
as deeply virtual ||, Q or wide-angle || [7j Compton scattering. The latter is in fact related 
to two-photon annihilation by crossing. It is important to realize that, since we take the BB 
system to have large invariant mass, the qq — ► BB transition can only be soft if the additional 
qq and possibly gluon pairs created in the hadronization process have soft momenta. In other 
words, the initial quark and antiquark must each take approximately the full momentum of 
one final state hadron. Compared to the case of mesons, our derivation for baryons will make 
the additional, plausible, assumption that the process is dominated by configurations where it 
is the quark that approximately moves in the direction of the baryon, whereas the antiquark 
approximately moves in the direction of the antibaryon. This assumption is equivalent to the 
valence quark approximation, widely used in other contexts. 




Figure 1: (a) Handbag factorization of baryon pair production 77 — > BB at large s, t, and u. 
The hard scattering subprocess is shown at leading order in a s , and the blob represents the 
BB distribution amplitude. A second graph is obtained by interchanging the photon vertices, 
(b) The physical mechanism of the handbag diagrams. The quark hadronizes into B and the 
antiquark into B, with any number of soft partons connecting the two parton-hadron vertices. 

For both processes, wide-angle Compton scattering off baryons and two-photon annihilation 
into BB pairs, the handbag contribution can dominate for large but not asymptotically large 
Mandelstam variables. Asymptotically the leading-twist contribution will take over, where in 
contrast to the handbag mechanism all valence quarks of the involved hadrons participate in 
the hard scattering ||. The handbag contribution formally represents a power correction to 
the leading-twist one. The onset of the leading-twist regime is however expected to occur for s 
much larger than experimentally available. A more detailed discussion of the relation between 
the leading-twist and soft handbag mechanisms, and of other power suppressed contributions 
is given in . 
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Two-photon annihilation into hadrons pairs has also been studied for the case where one of 
the photons has a virtuality Q 2 much larger than the squared invariant mass s of the hadron 
pair. In this kinematics, which is complementary to the one studied in the present article, 
handbag factorization of the process amplitude has been shown to hold for asymptotically 
large photon virtualities M. In other words, the handbag provides the leading-twist 

contribution in the limit of large Q 2 at fixed s. 

Our paper is organized as follows: In Sect. |2] we define the BB distribution amplitudes and 
discuss some of their properties. Sect. |3| is devoted to the calculation of the handbag amplitude 
for 77 — > BB. Flavor symmetry properties of the handbag amplitude are investigated in 
Sect. f|, and a comparison to experiment is presented in Sect. ||. We end our paper with a few 
concluding remarks. 

2 Distribution amplitudes for baryon-antibaryon pairs 

Generalized distribution amplitudes for meson pairs have been discussed in detail in the 
literature || |Hj. Here we introduce their counterparts for baryon-antibaryon pairs and 



present some of their general properties. We use light-cone coordinates v = [v + ,v , v_J with 
v ± = (v° ± f 3 )/y / 2 for any four- vector v, and define BB distribution amplitudes in light-cone 
gauge by 



%~ e~ izP+x ~ out (B(pu)B(p'u') | g(x) 7 + g(0) | ) = 

P+ _ 

®v( z > C, s) u(pv)^ + v{p'u') + $ q s (z, C, s) — u(pv) v(p'v') , 
^f- e~ lzP+x ~ ont (B(pv)B(p'i/) | g(a;)7 + 75 ?(0) | ) = 

P+ _ 

$a( z , C, s) w(pz/)7 + 75 v(p'u') + ® q P (z, C, s) — u{pv)^v{p'v') (1) 

with x = [0,£~,OjJ. Here m denotes the mass of the baryons and z/, v' their helicities. We 
have further introduced the sum P = p + p' of the baryon momenta, the invariant mass s = P 2 
of the baryon pair, and the skewness 

C = P + /P + - (2) 

In the following we will also use the notation ( = l — We have not displayed the dependence 
of the BB distribution amplitudes on the factorization scale /i 2 , which is governed by the well- 
known evolution equations for the distribution amplitudes of a single meson with appropriate 
quantum numbers. 

The BB distribution amplitudes are the time-like versions of generalized parton distribu- 
tions for a baryon B. Let us briefly comment on the relation of our definitions ([l|) with those 
of the generalized parton distributions H, E, H, E, introduced in M. Comparing the Lorentz 
structures that multiply the distributions and taking into account that p' + turns into —p' + 
under crossing of B(p'u') to S(pV), we recognize $,4 and $p as the respective counterparts 
of H and E. In the vector channel one may use the Gordon decomposition 

i 1,1 
u(pv)7: — a +p (p + ep')pv(p'u') = -(1 + e) u(pv)^ + vip'u') - - — (p - ep') + u(pv)v(p'v') (3) 
2m 2 2m 
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with e = ±1 to trade the scalar current for the tensor one. By crossing the defining relation for 
H and E one would obtain the scalar current u(pv) v{p'v') multiplied with {p'—p) + = (1—2()P + 
instead of P + . Defining a distribution amplitude with such a prefactor would however 
introduce an artificial singularity of $5(2, (, s) at £ = |, since there is no symmetry by which 
(p' — p) + $s(z, (, s) = P + $5(2, C) s ) has to vanish at p + = p' + . This is in contrast to the case 
of the generalized parton distribution E, where due to time reversal invariance the product 
{p' — p) + E occurring in its definition is zero for p + = p' + . 

Integrating ([D) over z reduces the bilocal BB matrix elements to local ones. In analogy to 
the space-like case we obtain a set of sum rules, 

F?{s) = Pdz&i&Cs) fori = V,A,P, 
Jo 

(i-20iioo = f'dzn^cs) (4) 

Jo 

with time-like form factors defined as 

ont (B(pv)B(pV)\q(0)Yq(0)\0) = F v u(pu)Y v(p'p') + Fj {p> ' u(pu) v(p'v') , (5) 

out ( B{pv)B{p'v') I g(0) 7 p 75 <?(0) | > = F% u{pv)l p ls v(p'u') + F q P W+Pr u(pu) l5 v(p'u') 

for each flavor. Appropriate combinations give the form factors of the electromagnetic and 
weak currents, for instance the magnetic and Pauli form factors 

Gm(s) = £ e q F v (s), F 2 (s) = £ e q F$(s). (6) 
<i q 

The relations (§) are valid for any physical value of the skewness (. They also hold for any 
value of the factorization scale fi 2 of the distribution amplitudes, since the vector and axial 
vector currents have zero anomalous dimension and the form factors Fy,s and Fa,p are scale 
independent. Taking higher moments in (2z — 1) leads to \j? dependent form factors of local 
operators with derivatives, multiplied by polynomials in (2( — 1). 

From charge-conjugation invariance we find the symmetry relations 



(z,(,s) = $?(*,C,s) iori = V,A,P, 

n( Z ,c,s) = -n(z,c,s) (?) 

with z — 1 — z. For £ = | the distribution amplitudes &v,a,p are hence symmetric under the 
replacement z z, while $5 is antisymmetric (but not zero). One may consider I?I?-states of 
definite charge-conjugation symmetry 

\C±) = \ \B(p, v) B(p'v') t B(p', v') B(pu)) , (8) 

satisfying 

C I C± ) = ± I C± ) , (9) 

where the operator C implements charge conjugation in Hilbert space. Replacing the state 
(B(p,v) B(p'u') I with (C±| in the definition ([I]), we obtain on its right-hand side the linear 
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combinations 



for i 
for i 



v,s, 

A, P. 



(10) 
(11) 



where we have used the symmetry relations (|7|). With the same relations one finds that $ 



and & 9 sap are °dd under the replacement C ^ C> which implies zeroes of these distribution 
amplitudes at £ = ~. Note also that the C-even combinations ^ys are antisymmetric under 
the replacement z <-> £ and therefore disappear in the sum rules (Q). This is consistent with 
the properties of the form factors Fy and Fg, which are C-odd. The reverse situation occurs 
for the C-odd combinations , which do not enter (|]) in agreement with the C-even nature 
of Fa and F P . 

The distribution amplitudes $j(z,C, s) are complex quantities, with phases due to the 
interactions in the BB system. Because of time reversal invariance they also parameterize 
matrix elements with baryons in the initial state, 

(h '~ izP+x ~ (0\q{0)j + q(x) | B{pv)B{p'v') ) in = 



«(+) 
V 



2tt 



P+ 



® v (zX,s)v(p'is')7 + u(pis) + C, s) — v(p'u') u{pv) 

2m 



dx 
~2T 



izP + x 



(o|g(o) 7 + 75g(^|5(H£(A')>in : 

p + 



z,(,s)v(pv')>y + >y 5 u(pv) 



2m 



v(p'u')j 5 u{pu). 



(12) 



Notice the change of sign in front of $p. 

We finally remark that distribution amplitudes for pairs B\B 2 involving different baryons 
can be defined as in ([[]). The charge conjugation relations (0) do not hold for these quan- 
tities, but they do for distribution amplitudes of the symmetrized states |i?i(pz/) £> 2 (pV) + 
B 2 (pu)B 1 (p'u')). 



3 The handbag amplitude 



We will now derive the expression for the soft handbag contribution to 77 — > BB. The first 
steps of the derivation go in complete analogy to the case of meson pair production. We thus 
start by summarizing the corresponding results of 0, and then proceed from the point where 
the different nature of baryons and mesons leads to important differences. In |l| we found an 
appropriate frame to be the cm. of the reaction, with axes chosen such that the process takes 
place in the 1-3 plane and the outgoing hadrons fly along the positive or negative 1-direction. 
Thus, we have baryon momenta 



V 



1,1, V2f3e 1 



P 



s r 



1,1, -V2(3e 1 



(13) 



with the relativistic velocity (3 
The photon momenta read 



q 



1 — 4m 2 / 's and ei = (1, 0). We hence have skewness £ 



- 1 + sin 6,1— sin 6 , \[2 cos 6 ei , 



2' 
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s r 



1 — sin 6,1 + sin 9 , — y/2 cos 9 ei 



(14) 



where 9 is the cm. scattering angle. In terms of the usual Mandelstam variables we have 



cos 9 



t — u 



sin 9 



2Vtu 



(15) 



up to corrections of order m?/s. The handbag amplitude for our process can be written in 
terms of the hard scattering kernel for 77 — > qq 



H af} {k,k') 



(k - q) • 7 , , {q-k')-j 
e -lT, no , • e "7 + e -77 , ao , ■ 6 '7 



(/c — q) 2 + ie 



(q — k') 2 + ze 



(16) 



a/3 



with photon polarization vectors e = e(q, fi) and e' = e(q', //), and a matrix element describing 
the qq — > BB transition. Our starting expression thus is 



A = £( ee «) 2 Jd'k mt (B(pv) B(p'u') I Tq a {x) q p (0) |0) H a0 (k, k>) 



(17) 



where the summation index q refers to the quark flavors u, d, s. In H we have omitted terms 
suppressed by the current quark masses. 

As discussed in detail in the qq — * BB transition at large invariant mass s can only be 
soft if the incoming quark and antiquark have small virtualities and each carry approximately 
the momentum of the baryon or antibaryon. To quantify this, we define z = k + /P + and 
parameterize the on-shell approximations of the quark and antiquark momenta k and kl as 



k 



s r 



k' 



z, z, -V2 



z, z, v2zz ej_ 

where e_i_ = (cosy?, sin <p). The requirements derived in [|lj then read 

A 2 

2z — 1, sin ip ~ — , 



zz e_|_ 



18) 



(19) 



where A is a hadronic scale of order 1 GeV. In addition, the minus- and transverse momenta 
oik — k and k' — k' must be of order A 2 / As shown in [0], 0], the dominant Dirac structure 
of the soft matrix element in (|i~7D involves the good components of the quark fields in the 
parlance of light-cone quantization. Projecting these out we have 



A , , 



d 4 k 



V4k+k'+ 



H^(k, k') S q (k, k') + HL,(k, k') S 5 q (k, k') 1 + 0(A 2 /s 



(20) 

where we have now made explicit the dependence on the photon and baryon helicities fi, // 
and u, z/, respectively. Here we have introduced the soft matrix elements^] 



d A x 



-ik-x 



out 



5(H5(A')|rg(x)7 + g(0)|0), 



(21) 



4 Note that we define soft matrix elements for states with definite baryon momentum here, and not for states 
(C+| with definite charge parity as in Q. 
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and S% with 7^ 



, replaced by 7 + 75. The hard subprocess amplitudes of 77 

sum and difference of the quark read 



qq for the helicity 



A=±l/2 

J2 2\u(k\)H^,(k,~k')v{k',-\) 

A=±l/2 



(22) 



where we have approximated the parton momenta with their on-shell values. The expressions 
(EO) and (|2"2] ) imply the phase conventions for light-cone spinors given in the appendix. With 
these conventions the behavior of helicity amplitudes under a parity transformation is A vv i^ = 



where 77 is the product of the intrinsic parities of the four particles 



involved. 

According to our hypothesis that the qq — > BB transition is dominated by soft processes, 
the matrix elements S q (k,k') and Sg(k,k') should be strongly peaked when (|T9|) is fulfilled. 
Depending on whether ip pa or ip pa tt this means that we have k pa p or k' pa p. The case 
k' ~ p corresponds to the antiquark hadronizing into -B and the quark into S. This requires 
sea quarks with very high momentum fraction in a baryon. We expect this to be disfavored 
compared with the case k pa p, both from phenomenological and theoretical considerations. A 
rather direct piece of information is for instance the ratio of quark and antiquark distributions 
in a nucleon at large momentum fraction x. Neglecting configurations with k' ~ p compared 
with k pa p, we proceed by Taylor expanding H(k, k ) and TC 5 (k, k ) around z — | and <p = 0, 

We get H<lI = 7if + 



keeping only the leading order in A 2 /s 
H+-(k,k') = W_ + (jfe,jfc / ) 

H 5 + _(k,k') = -H 5 _Jk,k') 



u/t — yi/u 
^u/t + Ji/u 



and 

„ * - 



+ 0(A 2 /s) 
+ 0(A 2 /s), 



(23) 



where t and u are the Mandelstam variables of the hard subprocess, with quark and antiquark 
momenta put on shell. For better legibility, explicit helicities are labeled only by their signs 
here and in the following. The amplitudes with equal photon helicities \i — // will be nonzero 
at next-to-leading order in a s , in analogy to the photon helicity flip transitions in large-angle 
Compton scattering fTI |. With (|23| ) the loop integration in ( p0|) now only concerns the soft 
matrix elements and leads to moments of BB distribution amplitudes. With dk + / \fAk + k' + ~ 
dz and the definitions of Sect. 0we obtain 



A 



dz \ H^(p,p') 



+ $ 5 



2' 



2 ■ ' 



' 2P 

1 

4m 
1 



u(pu)^f + v(p'v') 

u{pv) v{p'u') 
-M(pz/)7 + 7 5 v(pV) 



(24) 



is) 



2P 

— u(pu)j 5 v(p'u') 



+ 0(A 2 /s) 



Because of (0) the scalar distribution amplitude $5 decouples in our frame with ( 
ing the spinor products with the conventions given in the appendix, including terms suppressed 



|. Evaluat- 
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only by mj yfs, we arrive at our final result for the 77 — > BB amplitudes: 

s 




( 1) <A— V— V 

t — U 



Aira. 



elm 



R v (s) + 2v8 v - vl R A (s) + R P (s 



where we have defined the annihilation form factors by 



Ri(s) 



s) 



— Sw> Rp{s) 
2m 



for i — V,A,P . 



(25) 



0(A 2 /s) 



(26) 



with F? from (|j). As in wide-angle Compton scattering off baryons there are only three 
independent form factors. In Compton scattering it is the pseudoscalar rather than the scalar 
form factor that does not contribute, due to different choices of the reference frames. The 
unpolarized differential cross section is given by 



da 



BB) 



1 



s 2 sin 6 
Several comments on our result 



Rv(s] 



COS" 



9 + 



R A (s) + Rp(s) 



+ 



4m 2 



Rp(s) 



(27) 

are in order. Unlike R A and Rp, the vector form factor 
Rv projects on the C odd part of the BB state, whereas a 77 collision produces of course its 
C even projection. This is a result of the approximations that lead from ( ^0[) to (|24D , namely 
of neglecting configurations where the q emerging from the two-photon annihilation hadronizes 
into the baryon B instead of the antibaryon B. To take this contribution into account, one 
can split the loop integration over = (k 1 ,^) into the two hemispheres where k 1 is either 
positive or negative. In the former case one can expand the hard scattering amplitudes around 
(p = 0, and in the latter around if = ir. For ip = tt we then get TC+- ~ 2(u — t)/ytu instead 
of 2(t — u)j \rtu as in (p3|). The sum over both hemispheres thus gives a result proportional to 
'H+-{p ) p') times 

" dk~d 2 k ± sgn(A; 1 ) S q (k, k') , (28) 



instead of the integral 



j dk~d 2 k L Sg(k, k') 



(29) 



which we used to express our result in terms of the distribution amplitude $y. The integrated 
matrix elements (|28|) and (|29|) have opposite behavior under charge conjugation, but to the 



extent that the region k « p' gives a small contribution compared to the region k ~ p, their 
difference can be neglected. In our derivation we have preferred the form ( p9[ ) that leads to 
matrix elements of light-cone operators with well-known properties, at the price of a loss in 
accuracy which we do not expect to be critical. We remark that for baryon helicities v = v' 
both (|28D and (|29D are odd under z <-> ~z in our reference frame, which results in a zero 
contribution to the amplitude after integration over z. This can be shown by performing a 
charge conjugation followed by a rotation of 180° around the 3-axis. The amplitudes with 
opposite baryon helicities v = —v' do not decouple in this way, because charge conjugation 
exchanges the helicities of B and B. We finally emphasize that the amplitude ( p5|) for the 
77 — > BB process does have the correct behavior under charge conjugation, which for our 
spinor convention reads A vv ',^' = (— l) u ~ v A u i v y^ in this channel. 
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For the sake of comparison let us mention what happens if we make the corresponding 
approximations in wide-angle Compton scattering. Using that contributions where a fast 
antiquark is emitted from and reabsorbed by the baryon are small, one may count them with the 
"wrong" sign in the Compton form factors Ry{t) and Rr{t) of 0. Replacing then the explicit 
factors 1/x with 1, we obtain approximations Ry(t) ~ J2 q e l F?(t) an d Rrif) ~ J2 q e 2 F%{t) i n 
terms of the space-like Dirac and Pauli form factors, in analogy with our result here. 



The amplitude ( |2q ) shows important differences compared with the one we obtained in |l | 
for production of a pair of pseudoscalar mesons. Similarly to the BB amplitudes with v = u', 
the matrix element corresponding to (p8|) vanishes there when integrated over z. Due to parity 
invariance the corresponding contribution from is also zero, so that the leading term in 
the Taylor expansion ( ^3l) of the hard subprocess gives a vanishing scattering amplitude. We 
thus had to include the first nonleading term in this expansion, which is proportional to z — z. 
If one were to do the same in the BB case, one would get a further term in (p5|). It would 
involve / dz(2z — 1) $5(2, |, s), which is a form factor of the quark energy-momentum-tensor, 
in analogy to the meson pair case. Such a subleading contribution would behave as 1/ (tu) in 
the amplitude instead of 1/y/tu, and give a 1/ sin 4 9 rather than 1/ sin 2 9 dependence in the 
cross section. 

Returning to our result (|25|) , we observe that the pseudoscalar annihilation form factor 
Rp generates amplitudes with equal helicities of the baryon and antibaryon, i.e., with their 
spin projections along p coupled to zero. The spins of quark and antiquark in the qq — > BB 
transition do then not sum up to those of the BB system, which implies that parton configura- 
tions with non-zero orbital angular momentum along p are required in Rp. One expects that 
at large s such quantities are suppressed compared with R v or Ra- An analog in the space- 
like region are the form factors Rr{t) vs. Ry(t), and their electromagnetic counterparts F 2 (t) 
vs. Fi(t). Recent measurements from Jefferson Lab [|EJ indicate that for —t between 1 and 
5.6 GeV 2 the ratio F 2 (t)/Fi(t) approximately scales as m/\f—i. Assuming a similar behavior 
of Rp(s) / Ra{s) one finds that the term with s\Rp(s)\ 2 in (p7|) will not start do dominate over 
the other terms with increasing s. 

We see in ( p7|) that the form factor Ry can be separated from the two others through 
measurement of the angular distribution of the BB pair, given data of sufficient accuracy. 
This requires some lever arm in 9, but must stay within the validity of our approach, which is 
not applicable for 9 near or it, where the 77 — > qq subprocess is no longer hard. 

A separation of Ra and Rp can only be performed with suitable polarization measurements. 
Our amplitudes Q25D are evaluated for states with definite light-cone helicities, which is natural 
within our framework and leads to simple expressions. In the unpolarized cross section this 
does not matter, but for polarization observables the use of the ordinary c.m.s. helicity basis is 
more convenient. The transformation from one helicity basis to the other can be found in |TB 
In our kinematics, the c.m.s. helicity amplitudes M. read 



+ 



m 



2vA- 



+ 2v'A v - 



0(m 2 /s) . 



(30) 



An observable capable of separating Rp from the other form factors is, for instance, the helicity 
correlation Cll between baryon and antibaryon, given by 



C 



LL 



da(++) - da(+-) _ 


Ra(s) 


+ Rp(s) 


2 + cos 2 9 


Rv(s) 


2 

s 

Am? 


Rp(s) 


2 


da(++) + da(+-) 


Ra(s) 


+ R P (s) 


2 

+ cos 2 9 


Rv(s) 


2 

_| s 

4m 2 


R P (s) 


2 



(31) 
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where da(u, v') is the cross section for polarized BB production. 

One may also consider the time-reversed process BB — > 77, which for the case of proton- 
antiproton collisions may be experimentally accessible and has already been mentioned in [H . 
Time reversal invariance relates the amplitudes of both processes by 

•A f i[i',vu' 1 = ^w',^' • (32) 

Up to an extra (1 — 4m 2 /s)~ 1 in the phase space factor, BB — > 77 has therefore the same cross 
section as 77 — > BB. The relation between the BB distribution amplitudes for baryons in 
the initial or in the final state has already been given in Sect. ^]. As in the case of wide-angle 
Compton scattering |15||, one can finally generalize our approach to the case of virtual photons, 



provided their virtualities are at most of order s. 



4 Flavor symmetry 

We are now going to discuss various BB channels where B is a member of the lowest-lying 
octet of baryons. We shall derive relations among the various amplitudes and form factors 
in order to simplify the analysis of experimental data on these cross section, and to explore 
generic consequences of soft handbag dominance. 

Relations among the various BB channels are obtained by exploiting flavor symmetry, i.e. 
isospin and [/-spin invariance. The latter is the symmetry under the exchange d <-> s, and 
relates for instance the pp and the S + S~ channels. Since the photon behaves as a [/-spin 
singlet while (p, S + ) and (£~,p) are doublets, [/-spin conservation leads to 

A(E + E~) ~ A{pp) . (33) 

In contrast to isospin breaking, which is known to hold on the percent level and will be neglected 
here, [/-spin violations cannot numerically be ignored. This is indicated in (B3) and in later 



relations by the approximate symbol. In analogy to (|33"D one has 



A(~-E + ) ~ A{E~E + ) . (34) 

Other consequences of [/-spin symmetry hold for the [/-spin triplet (n, ±[£° + v^A],- ) and 
the [/-spin singlet | [A — v^3S°] . Together with the corresponding transformation properties of 
the antiparticles one obtains 

A(E°E°) ~ A(nn) ~ i U(S°S°) - 3^4(AA) 



A(AS ) ~ ^(S°A) ~ 1— [A{KK) - ^(E U E U ) j . (35) 

Notice that the preceding [/-spin relations hold in any dynamical approach respecting SU(3) 
flavor symmetry. 

To proceed, we use that the handbag mechanism involves intermediate qq states, 

77 — > qq — > BB , (36) 
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and generically decompose the 77 amplitudes as 



A(BB) = J2e 2 q Af. 



(37) 



We have omitted helicity labels since the following results hold for any set of helicities. The 
decomposition (j37|) already follows from ([H]) and thus is more general than our result (p5|). 
It does not rely on our neglect of various 0(A 2 /s) effects, nor of contributions where a fast q 
hadronizes into a baryon. 

A characteristic feature of the handbag mechanism is that the intermediate qq state can 
only be coupled to isospin 1 = and 1=1, but not to / = 2. This leads to particular strong 
restrictions in the £ sector, where it reduces the number of independent partial amplitudes 

to three, one for each flavor. The absence of the isospin-two component of EE implies the 
following relation for the amplitudes: 



(3* 



which provides bounds on the (integrated or differential) E°S° cross section 



a(E+£-) - Ja(E~E 



< vV(£°x°) < 2 ( vM s+ ^ 



cx(£-£- 



(39) 



This follows from isospin invariance alone and thus is a robust prediction of handbag domi- 
nance. 

Combining the relations due to isospin and to [/-spin, in particular the absence of final 
states with 1 = 2 for the handbag, we find that all BB channels are described by only three 
independent partial amplitudes, which one may take to be those of the pp channel, A p , A p d and 
A p s . The compiled relations of the other partial amplitudes to the proton ones read 



An 
n u 


= 


An _ Ap 
A d — A u > 




K = 


A p , 




= Af^ 


A p > = A 1 ^ ~ 




aT = 


A*~~ 




= Af K 


~[A P + A P ], 






-Af 


A A 


- A A ~ 

— n d — 


-\[K + ^ d + Ai], 




A A ~ 

a — 


-\[2A 


.AS" 
A u 


- -A A ^° 




J 


Af°- 


-- 0, 




— n d 






Af* = 


= 0, 


Af 


= A d - 


A p , Af = Af ~ 




Af = 


Af ^ 



A p . 



~A p d , 



(40) 



We now take recourse to valence quark dominance, which allows us to use the ampli- 
tudes fl25|) and the form factors F^ B (s). Valence quark dominance implies F^' p {s) = 0. With 
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this simplification the symmetry relations ([JO) hold for the form factors F^' B {s) as well, sepa- 
rately for i = V,A, P. We emphasize that in the context of the soft handbag amplitude, valence 
quark dominance does not assume that non- valence Fock states are unimportant, since any 
number of soft partons with appropriate quantum numbers can connect the two parton-hadron 
vertices in Fig. |I|b. Rather, we neglect contributions from sea quarks that carry almost all the 
momentum of a baryon, which should be a good approximation. 

With regard to the accuracy of the present data on 77 — > BB we simplify further by taking 
a single value p for the d/u ratio of all proton form factors F v , F^, Fp, 

Ff* = pF?* t i = V,A,P. (41) 

For the numerical analysis we will perform in Sect. [| this is not a severe restriction, since we 
find the annihilation cross sections dominated by the sum of the axial and the pseudoscalar 
form factor. The parameter p in (|4l|) is then essentially the one for the combination Fa + Fp. 
For further simplicity we will assume p to be real- valued and independent of s in our analysis. 
The ansatz Ff ,p (s) oc F^ ,p (s) parallels the behavior of fragmentation functions for d — > p and 
u — > p transitions. The d/u ratio in fragmentation is not well-known; a value of p = 1/2 is 



for instance chosen in fLE] . For time-like form factors one obtains the same value by analytic 
continuation to the point s = 0, where the Dirac form factors are Fi' p (0) = 1 and F"' p (0) = 2, 
if one makes the assumption that their ratio does not change significantly between s = and 
large time-like s. A value of p = 1 on the other hand is suggested by the fact that both u 
and d are valence quarks of the proton, and in order to produce the proton two quarks have 
to be created from the vacuum in both cases. Still different, the Lund Monte Carlo event 
generator |L7| provides a value of only 0.25 for leading protons (with z > 0.5). 

On the basis of this simple model for the soft physics input to the handbag approach, we 
can write the BB amplitudes as 

A{BB) = r B A{pp), (42) 

where r B = ^b{p)- The ratios of differential or integrated cross sections are then determined 
by r\. These factors read 

l + 4p 
4 + p 

1+p 1 5 + 2p 

3 1 + 2p 
~ 2 4+p ' 

V3 l-2p 



1 \ 



T AVI 



AS — ' S°A 



2 4 + p 



rwo ~ , TV- — . 43 

4+p " 4+p V ; 

We recall that these relations receive corrections from SU(3) flavor symmetry breaking. An 
investigation of the size and pattern of such corrections is beyond the scope of this work. We 
will see that neglecting them at the present stage is justified by the accuracy of available data. 
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5 Comparison with experiment 



A suitable and sufficiently precise set of data would allow for an experimental determination 
of the annihilation form factors, quite analogous to the measurement of electromagnetic form 
factors. As already mentioned, the angular distribution of the BB pair allows one to separate 
Ry from the combination of and Rp given in fl27D . The present data [ IS ]- [ 21 1 on unpolarized 
cross sections (we exclude low-energy data from our study) does not permit such detailed 
investigation. Moreover, most data is taken at energies which are rather low for the kinematic 
requirement of large s, t, u in the handbag approach. Below yfs ~ 3 GeV the dynamics may 
be dominated by resonances. 

It has long been known M that for asymptotically large s the process is amenable to a 
leading- twist QCD treatment, where the transition amplitude factorizes into a hard scattering 
amplitude for 77 — > qqqqqq and a single-baryon distribution amplitude for each baryon. As 
already mentioned in the introduction, the leading-twist result J22| is way below the experi- 
mental data. This holds in particular if the single-proton distribution amplitude is close to its 
asymptotic form under evolution, for which there is growing evidence now j33| . In view of this 
we consider that we make an acceptably small error in our present work by altogether neglect- 
ing the leading-twist contribution to the processes in question. We remark that on the other 
hand the diquark model, which is a variant of the leading-twist approach, provides reasonable 
fits to the data, at least for the pp channel Notice that both the leading-twist and the 
diquark approach give real-valued amplitudes in the 77 annihilation channel. In contrast, the 
handbag approach makes no generic prediction: the phase of the amplitude is determined by 
the phases of the annihilation form factors, which may or may not be small. 

The annihilation form factors and the BB distribution amplitudes can presently not be 
calculated from first principles in QCD. Contrary to generalized parton distributions, they 
do not admit a direct representation as overlaps of light-cone wave functions ]7|, ^] either. 
Progress in describing generalized distribution amplitudes within a Bethe-Salpeter approach 
has recently been reported 
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No model calculation is currently available for the annihilation 
form factors in the s range where we need them. We will therefore determine these form factors 
phenomenologically. 

Let us start with information from other processes. The E760 collaboration [27| has 
measured the magnetic form factor of the proton in the time-like region for s in the range 
8.9 -7- 13.0 GeV 2 . For the scaled form factor a value of s 2 |C^f I — 3 GeV 4 has been found. 
With (H), (|26"D, (|4lD this implies for the annihilation form factor 



2.4 GeV 4 , 3 GeV 4 , 5 GeV 4 for p = 0.25, 



0.5, 1, 



(44) 



if we neglect the non-valence contribution from F v . As we already discussed, the form factor 
Rp involves parton orbital angular momentum. For lack of better information we estimate the 
magnitude of Rp by assuming 



2m„ 



m 



t Fi(t) 



2m p Ff(t) 



0.37 



(45) 



for large s and t, where the numerical value is from the measurement [ 12 1 of F^t) / Ff{t) in 
the range — t — 1 -r- 5.6 GeV 2 . 
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We now turn to the two-photon annihilation data. Integrating the cross section (E7I) over 
cos 9 from — cos #0 to cos 8 , we get 



a(77 BB) 



4 ™elm 




Ry{s) 



+ 



R B A {s) 



R b p(s) 



+ 



4m 2 



R b p(s) 



cos 6*o 



1? 



When comparing with data we need the integrated cross section for cos #0 
choice of the experiments, 



(46) 

0.6, following the 



a (77 -> BB) = 181nb GeV s 



1 




4m? 



Rf\ 2 + 0.134 |i$| 2 



(47) 



We fit this to the data on 77 — > pp above s = 6.5 GeV 2 ~ (2.55 GeV) 2 , trying to avoid as 
much as possible the region where the process is markedly influenced by resonances. Such a fit 
determines the combination of form factors in the curly brackets of (f|7]). Neglecting the term 
with Ry we get 

s 2 R P . s (s) = (6.5 ± 0.5) GeV 4 , (48) 
with the fit shown in Fig. |2], where we have introduced the abbreviation 

/ \ V2 

1 r>B\2 s I p B|2 ) 

i A + K P \ + - — — \K 



R 



off 



\R B + R B P \ 2 + 



Am 2 B 



(49) 



With our estimates (fPj) the contribution of Ry to the cross section (^) is at most 8%. Taking 
it into account would thus reduce RP eS by at most 4%, which is below the error in ([48]). If we 
further use the estimate (|45| ) of Rp/R A at s = 6.5 GeV 2 , we obtain 



s 2 \R P a(s) 



(4.9 ± 0.4) GeV 4 + (8.0 ± 0.6) GeV 



(50) 



DP ~ 



\R P A + R P P \, with an accuracy between 



where the errors are due to those in the fit (|48| ) and the range to the uncertainty of the 
relative phase between Ra and Rp. Using the same input we get the approximate relation 

s) ~ s~~ 2 is compatible with experiment in the s range we 
are investigating, a somewhat different falloff is not excluded by present experimental data. 
The s dependence of our fitted annihilation form factor coincides with the one predicted by 

6 of the 



Although a behavior of R p s ( 



as well as the corresponding behavior da/dt ~ s~ 



dimensional counting rules |28 
cross section (|27| ) at fixed angle 9. We emphasize that this does not imply the dominance of 
leading- twist contributions. It is also possible that, in a way similar to wide-angle Compton 
scattering 0, HJ, dimensional counting rule behavior is mimicked by soft physics over a large 
yet finite range of s. From our calculation of the handbag diagrams it is clear that the form 
factors appearing in (^) are only the soft parts of the matrix elements Ri{s) (i = V, A, P) 
defined by @ and (|26|) . According to general power counting arguments, the soft parts of Ry, 
Ra and \J~s~Rp will decrease faster than s~ 2 for very large s. The soft handbag contribution 
to the cross section da/dt then falls off faster than s -6 at fixed 9, and the hard leading-twist 
contribution will eventually dominate. We remark that in the spacelike region one can use a 
model based on wave function overlap to evaluate the soft parts of the Compton form factors 
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Figure 2: The scaled annihilation form factor s R^ s (s) as extracted from the data of |20 



using (|47| ) with i?y set to zero. The dashed line represents a fit to the data above 6.5 GeV . 



Rv{t) and -Ra(^) 0- Their asymptotic behavior in this model is a decrease like t -4 and only 
sets in for —t of order 100 GeV 2 . 

We observe that the annihilation form factor ( fi8|) is of similar size as the time-like magnetic 
form factor of the proton. The situation is thus similar to the space-like region, where the Dirac 
form factor and the form factors for wide-angle Compton scattering off the proton also behave 
similarly and are of comparable magnitude J?], |15| . Recall that the Compton form factors Ra 
and Ry are given by moments in x of generalized parton distributions whose respective forward 
limits are the polarized and unpolarized quark densities Aq and q. If one assumes that at large 
x these generalized parton distributions have the same x-dependence as their forward limits, 
up to a common factor f(x,t) for both distributions, one obtains |-Ra(£)| < Rv(t) at large t 
as a consequence of the positivity bound |Ag| < q. For generalized distribution amplitudes 
there is no such constraint, and our estimates ( fH] ) and (|50| ) suggest that one may indeed have 
> |-Ry( s )| f° r th e annihilation form factors in the s-range of our fit. 

Using the result (^) for R p eS we can now discuss the cross sections for other BB channels. In 



view of the large uncertainties of the data |L9, |21[ we do not attempt to include effects of flavor 



symmetry breaking and directly use the relations ( f43|) to investigate the relative strength p 
between dd — > pp and uu —>■ pp transitions. In Fig. |3] we show the cross sections for two-photon 
annihilation into AA and S°S° pairs, with the bands corresponding to the range 

p = 0.25 ^0.75. (51) 

According to our discussion in Sect. £| such values are physically quite plausible. Values of 
p significantly different from (|51|) are not favored by the data. The estimate ( |51]) should 
of course be interpreted with due care, given theoretical uncertainties induced by the rather 
low s-values of the data (the respective production thresholds for AA and S°E° pairs are at 
y/s = 2.23 GeV and 2.39 GeV), the assumed s-dependence of _Rg ff in (0), and the simplifying 
assumptions that give the relations (f|3|) between BB channels in terms of a single real-valued 
parameter p. We recall that the value in ( |51| ) essentially refers to the form factor combination 
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Figure 3: The cross sections for two-photon annihilation into AA (left) and pairs (right) 
versus W = ^/s. Data is taken from |19|, |2l|. The bands correspond to the range p = 0.25-^-0.75 
in conjunction with the form factor B%. s (s) from (^), as explained in the text. 



_Rg ff ~ [R^ + Rp\, which dominates the integrated cross section (|47"D . Other BB channels can 
now easily be predicted from (|43|) . A special role is played by the mixed channels AS and 
£°A, whose cross sections vanish at p = 1/2. For these the estimate fl5l|) provides an upper 
bound 2 

<x(AS + S°A) ~ ~ *(pp) < 0-02 a{pp) , (52) 

whose precise value should again be taken with care, given our discussion above. 

In Fig. [| we show the angular distribution for 77 — > pp. Unfortunately, data exists only 
for rather small energies, where our kinematical requirements that —t and — u should be large 
compared to, say, the squared proton mass, can hardly be met. Furthermore, the influence of 
resonances may not yet be negligible, for which there might be a little hint in Fig. f|. At the 
energy of ^/s ~ 2.3 GeV (not shown in the figure) the data of JIB, ^] exhibit a maximum 
at 90°, which is a clear signal for the dominance of low partial waves and may be due to 
resonances. The comparison of the handbag result with the available data should therefore 
be interpreted with due caution. The curve in Fig. |4] shows the angular distribution of the 
handbag result when R v is neglected in (|27D. Taking our estimate ([H|) of R v for p = 0.5, 
together with the result ([48]) for _R^ ff , we get a change in the distribution that is too small to 
be seen in the figure. If on the other hand we take the value of Ry which corresponds to p = 1 
in (0), the angular distribution becomes somewhat steeper. We emphasize however that the 
region where —t or — u is smaller than 1.5 GeV 2 corresponds to | cos#| > 0.5 for the data in 
Fig. In this region the handbag result has to be taken with more than a grain of salt. 



6 Concluding remarks 

We have discussed the handbag contribution to two-photon annihilation into baryon-antibaryon 
pairs at large energy and large momentum transfer. Our main result is to write the amplitude 
as a product of a parton-level amplitude for 77 — > qq and annihilation form factors given by 
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Figure 4: Normalized angular distribution of 77 — > at s = 7.3 GeV 2 . The solid line 
represents the handbag result (^) with the form factor R v neglected. Data is taken from 



moments of the BB distribution amplitudes. In our derivation we have to explicitly neglect 
contributions where the antiquark nearly takes the momentum of the baryon and the quark the 
momentum of the antibaryon. On the other hand, quark off-shell effects in the hard scattering 
and the bad components of the corresponding field operators are shown to be suppressed 
parametrically. An alternative treatment of the processes under investigation is possible using 
double distributions P§[| . Our results also apply to the annihilation process pp — > 77, whose 
form factors and amplitudes are related to those for two-photon annihilation by time reversal. 

The factorization of the soft handbag diagrams is analogous to the one in wide-angle Comp- 
ton scattering. For the latter it has been shown that this factorization remains valid when 
taking into account next-to-leading corrections in a s to the parton- level subprocess [|ll|], and 
one may expect that the same holds for the time-like processes considered here. 

The handbag contribution formally represents a power correction to the leading-twist hard- 
scattering mechanism, but it seems to dominate at experimentally accessible energies. We find 
that the data for various BB channels is compatible with annihilation form factors approx- 
imately behaving as 1/s 2 for s between 6 and 12 GeV 2 , a counting rule behavior typical of 
many exclusive observables. Fitting the form factors to the data, we find that for protons 
the sum of the axial and pseudoscalar annihilation form factors Ra + Rp is dominant and 
somewhat larger than the time-like magnetic form factor. A further test of our approach is 
the approximate 1 / sin 2 9 angular dependence of the cross section, which agrees rather well 
with the VENUS data. According to our estimates, the Ry term with its additional cos 2 9 
dependences is likely too small to be seen in the presently available data. Flavor symmetry 
and the absence of / = 2 components in the qq intermediate states relate pp production to 
the BB channels where B is a member of the lowest lying baryon octet, up to presumably 
moderate effects of flavor SU(3) breaking. Fixing the relative strength p of the form factors 
governing dd — > pp and uu — > pp transitions from suitable and sufficiently accurate data of two 
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BB channels allows one to predict all other ones. 

We emphasize that our comparison with experiment suffers from the low energies where 
data is currently available. For these energies the kinematical requirements of the handbag 
approach are hardly satisfied. Nevertheless we arrive at a satisfactory description of the data 
for the three channels pp, A A and £°£°, taking as soft physics input the effective form factor 
i?g ff (s) and the flavor parameter p with values in agreement with the physical interpretation 
of these quantities. We finally remark that measurement of the process pp — > 77 with better 
statistics and at higher energies would likely be possible at the proposed HESR project at 
GSI IBOl. 



Note added 

After this work was finished, new data for the pp channel was published by OPAL f53]| . In 
the high-s range it agrees with the results from CLEO and VENUS we have referred to in 
Section |5l. 
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Appendix: Spinor conventions 



In our calculations we have used spinors for (anti)quarks and (anti)baryons that correspond 
to states with definite light-cone helicity |HTJ. In the usual Dirac representation they read 



u(p, +) 



v(p, +) 



/ p° + p 3 + m \ 



N" 1 



p 1 + ip 2 
p° + p 3 
\ p 1 + ip 2 



m 



J 



N' 1 



P 



ip 



\ 



p" — p + m 
—p 1 + ip 2 
V p° + p 3 + m J 



u(p, 



v(p, -) 



N- 1 



N~ 



I —p 1 + ip 2 \ 

p° + p 3 + m 

1 • 2 
p — ip 

\ —p° — p 3 + m j 

I —p° — p 3 + m \ 

1 • 2 
— p — ip 

n s 
— p — p — m 

\ —p 1 — ip 2 



(53) 



where N = J2(p° + p 3 ). This corresponds to the phase conventions used by Brodsky and 



Lepage, cf . |32j , and also to those of Kogut and Soper [31| if one takes into account that they 
use a different representation of the Dirac matrices. The antiquark spinors in (|53D satisfy the 
charge conjugation relations v(p,u) = S(C)u T (p,u) with S(C) = ry 2 7°. For massless spinors 
one simply has v (p, v) = —u(p, —v). 
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